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1.  Introduction 


We  considered  the  following  scenario:  A  steel  projectile  perforates  a  sheet  of  titanium  armor 
and  is  heated  by  plastic  work  during  the  penetration  process.  The  hot  steel  residual  then  comes 
to  rest  in  a  bed  of  M30A1  propellant.  The  subsequent  temperature  rise  in  the  propellant,  if  of 
sufficient  magnitude  and  duration,  can  lead  to  initiation  (7). 

Our  goal  was  to  develop  a  computational  methodology  to  study  this  scenario.  Section  2 
describes  our  application  of  the  CTH  Eulerian  hydrocode  (2)  to  the  steel-titanium  penetration 
problem.  Section  3  describes  our  application  of  EnSight  (3)  to  the  extraction  from  the  resulting 
CTH  Spyplot  output  files  the  geometry  of  and  temperature  field  within  the  steel  residual. 

Section  4  then  describes  our  use  of  LS-DYNA  ( 4 )  to  study  heat  conduction  from  the  residual  to  a 
single  grain  of  M30A1  propellant.  Section  5  concludes  with  a  summary  and  assessment. 


2.  Titanium  Perforation  Simulation  With  CTH 

2.1  Problem  Definition 

A  right  circular  cylinder  (RCC)  of  4340  steel  impacts  a  titanium  plate  at  normal  incidence  and  a 
speed  of  433  m/s  (figure  1).  The  RCC  has  a  length  and  diameter  of  20  mm.  The  target  plate  has 
a  thickness  of  6.35  mm  and  a  square  face  with  127-mm  edge  length. 

2.2  Meshing 

The  stationary  (Eulerian)  mesh  consisted  of  6.24M  cells,  each  a  cube  of  0.635-mm  edge  length. 
In  terms  of  the  coordinate  system  defined  in  figure  1,  the  x-length  of  the  mesh  was  95.25  mm, 
and  the  y-  and  z-lengths  were  each  129.54  mm.  Symmetry  boundary  conditions  were  not 
employed. 

2.3  Material  Properties 

The  projectile  was  composed  of  martensitic  4340  steel.  The  target  plate  was  composed  of  the 
titanium  alloy  Ti-6A1-4V.  The  material  models  used  are  described  next.  The  values  assigned  to 
all  material  constants  are  given  in  table  1 .  For  most  parameters,  the  source  of  the  value  is  also 
given  in  the  table.  For  the  others,  there  is  more  uncertainty  in  the  value. 
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Figure  1.  An  RCC  of  4340  steel  strikes  the  center  of  a  square  titanium  plate  at  normal  incidence  and 
433  m/s.  (The  drawing  is  not  to  scale.) 

2.3.1  Equation  of  State 

The  dilatational  response  of  the  4340  steel  and  the  titanium  were  both  modeled  by  a  Mie- 
Griineisen  equation  of  state  with  a  constant  Griineisen  coefficient,  yo.  The  shock  Hugoniot  curve 
was  described  by  the  linear  relationship  between  the  shock  speed,  Us,  and  the  particle  speed,  tip, 

Us  =  Cs  +  Sup .  (1) 


Here,  Cs  is  the  sound  speed  and  S  is  a  material  constant.  Note  in  table  1  that  the  titanium  alloy 
was  assigned  properties  obtained  for  aluminum. 

2.3.2  Strength 

The  titanium’s  strength  was  represented  with  the  FCC  form  of  the  Zerilli- Armstrong  model  (5), 
according  to  which 


Y  =  C0  +  C,  exp(-  C36>  +  C46»ln^)+  C5en . 


(2) 
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Table  1.  Material  properties  in  the  section  2  CTH  calculations. 


4340  Steel 

Ti-6A1-4V 

p  (kg/m3) 

7850  (6,  7,  8,  9 ) 

2660  (5)  for  A1 

To 

2.17 

1.97  (5)  for  A1 

Cs  (m/s) 

4570  (< 8 ) 

5220  (8,  9)  for  A1 

5 

1.49  (9) 

1.37  (8,  9)forAl 

Cp  (J/kg-K) 

475  (6,  7) 

896  (<5)  for  A1 

V 

0.29  (7) 

0.30 

Co  (GPa) 

— 

1.116 

Ci  (GPa) 

— 

0.24  (10) 

C3  (1/K) 

— 

0.0024  (70) 

C4  (1/K) 

— 

0.0043  (10) 

C5  (GPa) 

— 

0.656  (10) 

n 

— 

0.50  (10) 

A  (GPa) 

0.72  (77) 

— 

B  (GPa) 

0.78  (77) 

— 

C 

0.004  (77) 

— 

M 

1.00  (77) 

— 

N 

0.106  (77) 

— 

0m  (K) 

1782  (6,  7) 

— 

&r  (K> 

294 

— 

D\ 

See  (72) 

See  (72) 

d2 

See  (72) 

See  (72) 

d3 

See  (72) 

See  (72) 

D 4 

0.002  (7) 

See  (12) 

d5 

0.61  (7) 

See  (12) 

Tfrac  (GPa) 

-2.5 

-5.9 

Here,  Y  is  the  yield  stress,  0  the  temperature,  ^the  equivalent  plastic  strain,  and  s  the  time  rate 
of  change  of  equivalent  plastic  strain  nondimensionalized  by  a  reference  rate  of  1/s.  The  model 
introduces  six  material  constants:  Co,  C\,  C3,  C4,  C5,  and  n. 

The  steel’s  strength  was  represented  with  the  Johnson-Cook  strength  model  (6),  according  to 
which 

Y  =  [A  +  B(sp)N\\  +  C\ns1^-(e*)M  \.  (3) 


In  equation  3,  6*  is  the  homologous  temperature,  defined  by 
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0-0r 
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(4) 
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where  Or  is  the  room  temperature  and  0m  is  the  material’s  melting  temperature.  The  seven 
material  constants  are  A,  B,  C ,  N,  M,  9m ,  and  0r .  In  addition  Poisson’s  ratio,  v,  is  used  to 
compute  elastic  strains. 

2.3.3  Damage 

The  Johnson-Cook  fracture  model  (7)  was  applied  to  both  4340  steel  and  Ti-6A1-4V.  This  model 
introduces  a  fracture  strain,  £/,  which  is  related  to  the  local  stress,  strain  rate,  and  temperature 
according  to 

ef  =  [.D,  +  D2  exp(z>3cr* )]  [l  +  Z)4  In e] [l  +  DS6*],  (5) 


where 


* 

C 7 


<J 


eq 


(6) 


Here,  P  is  pressure,  <req  is  the  von  Mises  stress,  and  s  and  s  are  again  equivalent  plastic  strain 
and  the  dimensionless  rate  of  change  of  equivalent  plastic  strain,  respectively.  The  Johnson- 
Cook  fracture  model  introduces  five  material  constants,  D\,  D2,  D3,  D4,  and  D5. 

In  addition,  CTH  offers  the  ability  to  insert  void  in  cells  that  exhibit  pressure  less  than  or  equal  to 
a  prescribed  value,  Pfrac.  This  feature  was  applied  to  both  metals,  with  the  P& ac  values  in  table  1 . 
The  negative  values  indicate  tension. 

2.4  Results 

Heat  conduction  was  neglected  in  the  CTH  perforation  calculations.  The  temperature  rise  in  a 
cell  during  a  given  time  step  was  calculated  by  dividing  work  increment  by  specific  heat.  The 
results  for  volume  fraction  and  temperature  of  steel  in  each  cell  were  written  to  Spyplot  files 
called  “spcth.”  These  are  the  solution  variables  needed  for  the  thermal  postprocessing  with 
LS-DYNA. 


Figures  2  and  3  show  the  regions  of  steel  and  titanium  at  intervals  of  20  ps.  These  regions  were 
identified  based  on  a  material  volume  fraction  of  0.5  or  greater.  The  problem  was  terminated  at 
200  ps  after  impact. 

2.5  Discussion 

The  problem  in  figure  1  was  studied  experimentally  by  Kyle  Bates  of  the  U.S.  Army  Research 
Laboratory’s  Survivability/Lethality  Analysis  Directorate.  In  the  experiment,  the  recovered 
titanium  plug  and  steel  residual  were  bounded  by  shiny  surfaces  (13)  indicative  of  plastic  flow 
and  possibly  adiabatic  shear  banding.  In  the  CTH  simulation,  damage  in  the  two  metals  was 
modeled  with  Johnson-Cook  fracture.  It  has  not  been  established  that  this  damage  model 
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Figure  2.  CTH  output  for  the  steel  vs.  titanium  at  t  =  0  to  100  ps  after  impact. 
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Figure  3.  CTH  output  for  the  steel  vs.  titanium  at  t  =  120  to  200  ps  after  impact. 
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introduces  effects  of  adiabatic  shear  localization.  An  effect  of  plug-bounding  shear  localization 
in  the  titanium  target  is  to  reduce  the  work  needed  to  perforate  the  target.  Thus,  the  computed 
plastic  work  and  temperature  rise  in  the  steel  projectile  are  likely  to  overestimate  the  actual  values. 

As  noted  in  section  2.3.1,  we  applied  to  the  titanium  alloy  values  for  equation  of  state  constants 
that  were  obtained  for  aluminum.  This  was  another  source  of  error. 


3.  Extraction  of  Residual  Geometry  and  Temperature  Field  Using  EnSight 


The  spcth  fdes  were  read  with  EnSight  (figure  4).  This  software  enabled  identification  of  those 
CTH  cells  that  exhibited  a  steel  volume  fraction  of  0.5  or  greater.  This  criterion  was  only 
satisfied  by  the  single  closed  volume  of  cells  shown  in  figure  5.  We  identified  this  region  as  the 
steel  residual. 

CTH  specifies  temperature  at  the  centroid  of  each  cell.  EnSight  was  used  to  interpolate 
temperatures  to  the  eight  vertex  nodes  of  each  cell.  The  temperature  field  is  shown  in  figures  6 
and  7.  The  former  shows  temperature  on  the  outer  boundary  of  the  steel  residual,  the  latter 
shows  temperature  in  the  interior  exposed  by  a  diametral  cut. 

EnSight  was  used  to  create  three  ascii  files: 

1 .  Element  cards 

Each  line  consists  of  nine  integers:  a  unique  number  assigned  to  each  CTH  cell  that 
satisfied  the  0.5  volume  fraction  steel  criterion,  followed  by  eight  numbers  identifying  each 
of  the  cell’s  vertex  nodes. 

2.  Node  cards 

Each  line  contains  four  quantities:  the  integral  node  number  and  the  x,  y,  and  z  coordinates 
of  that  node. 

3.  Initial  nodal  temperatures 

Each  line  consists  of  the  integer  identifying  a  vertex  node  deemed  to  be  part  of  the  residual, 
followed  by  its  initial  temperature  in  Kelvin. 

These  three  sets  of  cards  were  brought  into  LS-DYNA  Keyword  format  and  then  combined  to 
form  part  of  a  Keyword  input  file.  The  latter  was  then  brought  into  HyperMesh  (14),  a  software 
product  used  to  preprocess  for  LS-DYNA.  The  remainder  of  the  LS-DYNA  input  file  was  then 
created  within  HyperMesh.  This  remainder  consists  specifically  of  the  mesh  for  the  M30A1 
propellant,  the  descriptions  of  the  steel  and  M30A1  thermal  material  properties,  the  thermal 
boundary  conditions,  and  the  termination  time,  time  step,  and  thermal  solver  specification. 
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Figure  4.  Overview  of  computations. 


Figure  5.  The  CTH  cells  that  at  t  =  200  |as  are  identified  with  the  steel  residual. 
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Figure  6.  CTH  results  for  temperature  contours  in  the  steel  residual  at  t  =  200  ps;  these  temperature 
results  are  applied  as  initial  condition  in  the  LS-DYNA  analysis.  (The  fringe  levels  are  in 
Kelvin.) 


4.  Residual  and  Propellant  Heat  Conduction  Analysis  With  LS-DYNA 


4.1  Problem  Definition 

Horst  et  al.  (15)  describes  the  geometry  of  individual  grains  of  M30A1  as  RCCs  with  diameters 
of  10.6  mm  and  lengths  of  24.4  mm.  We  considered  a  single  grain  composed  of  a  homogeneous 
material  and  placed  in  direct  contact  with  a  hot  section  of  the  steel  residual.  A  Lagrangian  finite- 
element  (FE)  model  was  created  for  this  grain.  We  assumed  that  all  mechanical  deformation, 
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Fringe  Levels 
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Figure  7.  The  steel  residual  in  figure  6  is  sectioned  along  an  approximate  diameter  to  reveal  the 
internal  temperature  field.  (The  fringe  levels  are  in  Kelvin.) 

including  thermal  expansion  and  remaining  stress  waves,  are  negligible  in  the  subsequent  heat 
transfer  to  the  M30A1  grain  once  the  steel  residual  has  come  to  rest.  This  allowed  for  the  use  of 
LS-DYNA  in  an  implicit  heat  conduction  mode.  We  assumed  perfect  thermal  contact  between 
the  steel  residual  and  the  M30A1  grain. 

Let  Q,  and  Q2  be  the  regions  of  space  occupied  by  the  steel  residual  and  the  M30A1  grain, 
respectively.  Let  9X  (x,  t)  and  02  (x,  t)  be  the  temperature  fields  within  the  steel  residual  and  the 
M30A1  grain,  respectively.  Heat  flow  is  then  governed  by  the  equations 

— -  =  cqV2#! ;  x  e  Qj ,  t  >  0  (7) 

dt 

and 

dG~, 

- =  a2'V292 ;  x  e  Q2 ,  t  >  0,  (8) 

dt 
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in  which  a\  and  a2  are  the  thermal  diffusivities  of  steel  and  M30A1,  respectively.  These  are 
defined  by 


and 


ax 


P\cvi 


a2 


(9) 


(10) 


Here,  K\  ad  k2  are  the  thermal  conductivities,  p\  and  fh  are  the  densities,  and  cp i  and  cp2  are  the 
constant-pressure  specific  heats  for  the  steel  and  the  M30A1  grain,  respectively.  These 
conductivities,  densities,  and  specific  heats  were  assumed  to  be  material  constants,  independent 
of  temperature. 


The  initial  conditions  are 


6»! (x,0)  =  6»0 (x) ;  xeQ,  (11) 

and 

d2 (x,0)  =  #r  =  294K ;  x  e  Q2 ,  (12) 

where  <90  (x)  is  the  initial  temperature  field  within  the  steel  residual.  This  temperature  field  is 
identified  with  the  CTH  output  shown  in  figures  6  and  7.  6t  is  room  temperature. 


The  imposed  boundary  conditions  are  indicated  schematically  in  figure  8.  Along  the  interface 
between  steel  and  propellant,  perfect  thermal  contact  is  assumed,  so  that  temperature  and  heat 
flux  are  both  continuous. 


6X  (x,  t)  =  02  (x,  t) ;  x  e  Qj  u  Q2 ,  t  >  0 . 


K, 


d6>i(x,t) 

dnx 


=  -K, 


d02(x,t ) 
dn. 


;  x  e  Q,  uQ, ,  t  >  0. 


(13) 

(14) 


Here,  n\  and  n2  indicate  the  directions  of  the  outer  normal  to  the  boundary  of  the  steel  residual 
and  propellant  grain,  respectively.  Throughout  the  portions  of  the  steel  and  propellant 
boundaries  not  in  mutual  contact,  an  insulation  boundary  condition  is  assumed.  Let  ins  and 

SQ2  ins  be  these  insulated  sections  of  the  steel  and  propellant  boundaries,  respectively.  Then, 

=  0;  XGgQ  ;>q  (15) 

onx 

and 

56>2(x,0  =  Q.  xegQ  t>0  (16) 

on2 
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Figure  8.  Schematic  of  the  boundary  conditions  applied  in  the  heat  conduction  analysis  of  section  4. 
Temperature  and  heat  flux  are  continuous  along  the  residual-propellant  interface;  heat 
flux  is  zero  elsewhere  along  the  residual  and  propellant  boundaries. 

4.2  Material  Properties 

The  applicable  material  properties  for  this  heat  conduction  analysis  are  density,  specific  heat, 
thermal  conductivity,  and  thermal  diffusivity.  Only  three  are  independent,  with  the  fourth  related 
by  equation  9  or  10. 

Values  for  these  material  constants  in  the  case  of  4340  steel  are  listed  in  table  2.  The  source  of 
each  value  is  indicated. 


Table  2.  Material  constants  for  steel  and  propellant. 


4340  Steel 

M30A1 

p  (kg/m3) 

7823  (6,  7, 8,  9) 

2000 

cP  (J/kg-K) 

477.8  (6,  7) 

1100  (equation  10) 

v(W/m-K) 

38.11  (7) 

0.418(7(5) 

a  (m2/s) 

1.020  x  10“5  (equation  9) 

1.9  x  10-7  (7<5) 
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Miller  ( 16)  measured  thermal  conductivity  and  thermal  diffusivity  of  a  grain  of  M30A1  as 
functions  of  temperature  and  direction.  His  data  have  been  converted  to  SI  units  and  reproduced 
in  table  3.  “M30//”  and  “M30j_”  indicate  properties  along  and  orthogonal  to  the  axis  of  the  grain, 
respectively.  For  both  k  and  «,  measurements  were  generally  larger  along  the  grain  axis.  Since 
/rand  a  were  both  modeled  as  independent  of  temperature  and  direction  in  the  LS-DYNA 
simulations,  a  representative  value  for  each  quantity  was  used  (see  table  2).  Specific  heat  was 
then  determined  from  equation  10. 


Table  3.  Thermal  conductivity  and  diffusivity  of  an  M30  grain  (reproduced  from  Miller  [16]). 


M30± 

M30„ 

e  (K) 

Ar(W/m-K) 

a  (x  10  7  m2/s) 

v(W/m-K) 

a  (x  10  7  m2/s) 

254 

n.d. 

n.d. 

0.419  ±0.021 

2.34  ±0.09 

255 

0.329  ±  0.007 

1.78  ±0.03 

n.d. 

n.d. 

275 

0.352  ±0.010 

1.74  ±0.04 

n.d. 

n.d. 

276 

n.d. 

n.d. 

0.486  ±  0.029 

2.41  ±0.14 

295 

n.d. 

n.d. 

0.427  ±  0.025 

1.90  ±0.09 

296 

n.d. 

n.d. 

0.448  ±  0.025 

2.09  ±0.11 

297 

0.320  ±0.033 

1.46  ±  0.1 1 

n.d. 

n.d. 

321 

0.325  ±  0.022 

1.39  ±0.08 

n.d. 

n.d. 

322 

n.d. 

n.d. 

0.412  ±0.016 

1.88  ±0.04 

Note:  n.d.  =  no  data. 


4.3  Meshing 

The  Eulerian  mesh  of  the  CTH  output  was  converted  directly  into  Lagrangian  FEs.  Each  steel 
FE  was  a  cubic  8-node  brick  with  a  0.635-mm  edge  length  (figure  5). 

The  new  mesh  created  for  the  cylindrical  M30A1  grain  also  consisted  entirely  of  cubic  8-node 
bricks  with  a  0.635-mm  edge  length  (figures  9  and  10).  HyperMesh  was  used  in  creating  this 
mesh. 

An  end  face  of  the  M30A1  mesh  was  placed  in  contact  with  a  flat  region  on  the  periphery  of  the 
steel  residual  (see  figures  10,  11,  and  12).  This  flat  region  of  the  residual  was  comprised  of  a 
face  of  each  of  83  elements  and  formed  a  substantial  contact  surface  for  the  propellant  grain. 
Note,  however,  that  this  flat  region  of  the  residual  did  not  coincide  with  the  region  in  which  CTH 
determined  temperatures  to  be  the  greatest.  The  vertex  nodes  for  each  of  these  83  steel  element 
faces  shown  in  figure  1 1  were  shared  by  the  83  propellant  element  faces  shown  in  figure  10.  In 
this  way,  the  continuity  of  temperature  and  heat  flux  across  the  residual-grain  interface 
(equations  13  and  14)  was  imposed. 

4.4  Results  and  Discussion 

The  computed  temperature  fields  at  1  and  2  s  after  contact  are  shown  in  figures  13  and  14.  The 
latter  figure  is  a  diametral  section  to  reveal  the  internal  field.  Since  the  thermal  diffusivity  of 
steel  is  two  orders  of  magnitude  larger  than  that  of  M30,  heat  flows  more  rapidly  within  the  steel 
than  from  the  steel  to  the  propellant. 
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Figure  10.  The  entire  mesh  for  an  M30A1  grain.  A  face  of  each  of  83  M30  elements,  colored  red,  was 
contiguous  with  the  steel-residual  mesh. 
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Figure  11.  The  face  of  each  of  83  residual  elements,  colored  red,  was  contiguous  with  the  M30 
mesh. 
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Figure  12.  An  FE  model  for  a  cylindrical  grain  of  M30A1  was  generated  and  placed  in  contact  with  the 
steel  residual. 
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Figure  13.  Temperature  contours  on  the  outer  surfaces  of  residual  and  propellant  initially  and  at  1 
and  2  s  after  initial  contact.  (The  fringe  levels  are  in  Kelvin.) 
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Figure  14.  Temperature  contours  in  the  midplane  of  residual  and  propellant  initially  and  at  1  and  2  s 
after  initial  contact.  (The  fringe  levels  are  in  Kelvin.) 
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Figure  15  identifies  five  nodes  within  the  diametral  plane  and  at  the  elevation  (x  location)  of  the 
hottest  region.  The  separation  between  adjacent  nodes  is  the  element  size,  0.635  mm.  The 
temperature  histories  of  these  nodes  are  shown  in  figure  16.  Nodes  a  and  b,  in  the  interior  of  the 
steel  residual,  and  node  c,  at  the  residual-grain  interface,  cool  down  during  the  first  0.5  s  to  a 
uniform  temperature.  At  2.4  s,  when  the  calculation  was  stopped,  these  nodes  are  continuing  to 
gradually  cool  further.  Nodes  d  and  e  lie  within  the  propellant  grain  and  heat  up  over  time. 
There  is  more  disparity  between  the  temperatures  of  d  and  e  than  between  the  temperatures  of 
nodes  within  the  steel,  indicative  of  the  steel’s  greater  thermal  diffusivity. 


Figure  15.  The  five  nodes  at  which  temperature  history  is  provided  in  figure  16. 
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Figure  17  shows  the  initial  temperature  field  on  face  A,  the  face  of  the  propellant  grain  in  contact 
with  the  steel  residual.  Face  B  is  the  plane  of  propellant  material  recessed  one  element 
dimension,  0.635  mm,  from  face  A.  Nodes  1 A-6A  are  the  six  hottest  nodes  in  face  A  and  are 
shared  by  the  steel  residual.  Nodes  1B-6B  are  their  counterparts  with  the  same  x  and  z 
coordinates  but  on  face  B. 

Figure  18  shows  the  computed  temperature  histories  of  nodes  1A-6A.  These  nodes  begin  at 
similar  temperatures  and  cool  at  similar  rates.  Figure  19  shows  temperature  histories  for  nodes 
1B-6B.  Nodes  1B-6B  all  start  at  0r  and  subsequently  heat  up. 
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Figure  17.  The  initial  temperature  field  in  the  M30A1  grain,  and  the  six  nodes  at  which  temperature  histories 
are  provided  in  figure  18.  (The  fringe  levels  are  in  Kelvin.) 
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Figure  18.  Temperature  vs.  time  at  the  six  nodes  identified  in  figure  17  on  face  A  (the  surface  in  perfect 
contact  with  the  steel  residual),  and  ignition  times  reported  in  Rocchio  and  Wires  (7).  (Note 
the  logarithmic  temperature  scale.) 
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Figure  19.  Temperature  vs.  time  at  the  six  nodes  identified  in  figure  17  on  face  B  (the  M30  internal  plane 
recessed  0.625  mm  from  the  steel  residual). 


Figure  20  shows  temperature  histories  for  nodes  6A  and  6B.  Node  6A  begins  at  524  K.  At  2.4  s, 
temperatures  of  nodes  6A  and  6B  have  converged  to  328  and  308  K,  respectively. 

Figures  18  and  20  contain  data  points  labeled  “R&W  Lower”  and  “R&W  Upper.”  These  refer  to 
M30  ignition  data  reported  in  Rocchio  and  Wires  (7).  Rocchio  and  Wires  performed  temperature 
controlled  experiments  on  single  grains  of  M30.  For  a  range  of  temperatures  between  483  and 
533  K,  the  time  required  for  ignition  was  measured.  Their  results,  after  converting  from  Celsius 
to  Kelvin,  are  reproduced  in  table  4.  To  quote  from  Rocchio  and  Wires  (7),  “the  high  limit  was 
that  time  above  which  one  could  always  expect  sustained  combustion  and  the  low  limit  was  that 
time  below  which  sustained  combustion  did  not  occur”  (p.  14).  These  data  are  the  basis  for  an 
initiation  criterion  to  apply  to  the  LS-DYNA  temperature  solution.  In  figure  18,  since  all  points 
on  curves  6A  through  6B  lie  below  the  Rocchio  and  Wires  data  points,  one  can  conclude  that 
the  M30A1  grain  that  contacts  the  steel  residual  in  the  location  and  configuration  depicted  in 
figures  10-12  will  not  undergo  sustained  combustion.  Note  again  that  in  our  analysis  the  grain 
did  not  contact  the  residual  at  the  latter’s  hottest  location. 
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Figure  20.  Temperature  vs.  time  at  two  nodes  identified  in  figure  17,  one  on  face  A  and  one  on  face  B, 
and  ignition  times  reported  in  Rocchio  and  Wires  (7). 

Table  4.  Measured  times  required  for  sustained 
combustion  as  functions  of  temperature 
(reproduced  from  Rocchio  and  Wires  [7]). 


0(K) 

Slower  (s) 

Aipper  (s) 

533 

0.260 

0.378 

523 

0.574 

0.621 

513 

1.538 

2.329 

503 

3.422 

4.111 

493 

6.257 

7.045 

483 

9.400 

18.256 

The  issue  of  mesh  sensitivity  is  explored  in  the  appendix.  There  we  study  one-dimensional  heat 
conduction  between  a  semi-infinite  region  of  steel,  initially  at  a  uniform  non-zero  temperature,  in 
perfect  thermal  contact  with  a  semi-infinite  region  of  propellant,  initially  at  zero  temperature. 
First,  an  analytical  solution  is  obtained.  Next,  the  problem  is  simulated  with  LS-DYNA  using 
the  same  0.635-mm  cubic  elements  as  used  in  the  present  section  4.  Good  agreement  between 
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the  analytical  and  the  LS-DYNA  solutions  in  the  propellant  give  some  support  to  the  adequacy 
of  the  meshing  employed  here  in  section  4.  Note,  however,  that  the  problem  solved  in  the 
appendix  involves  a  uniform  initial  temperature  throughout  the  steel,  whereas  the  spatially 
varying  temperature  field  predicted  by  CTH  introduces  additional  length  scales  that  can  affect 
mesh  adequacy. 


5.  Concluding  Remarks 


5.1  Summary 

A  procedure  was  developed  for  performing  a  postprocessing  heat  conduction  analysis  from  a 
CTH  solution  for  a  residual  penetrator’s  shape  and  temperature  field.  LS-DYNA  was  employed 
for  the  heat  conduction  analysis.  The  interface  between  CTH  and  LS-DYNA  involved  the  codes 
EnSight  and  HyperMesh.  Figure  4  schematizes  the  overview. 

The  procedure  was  applied  to  perforation  of  a  titanium  plate  by  a  steel  RCC  at  433  m/s,  followed 
by  cookoff  of  an  M30A1  grain  by  contact  with  the  hot  steel  residual.  Heat  conduction  was 
neglected  in  the  steel-titanium  perforation  analysis,  performed  with  CTH.  Deformation  and 
motion  were  neglected  in  the  postprocessing  heat  conduction  analysis,  performed  with  LS-DYNA. 
The  temperature  field  as  a  function  of  time  within  the  M30A1  grain  was  compared  with  data  in 
Rocchio  and  Wires  (1)  for  initiation  time  as  a  function  of  temperature. 

5.2  Assessment 

The  heat  conduction  analysis  with  LS-DYNA  was  straightforward  and  presented  no  special 
challenges.  The  main  uncertainties  were  introduced  by  the  perforation  analysis  with  CTH.  In 
particular,  the  steel  and  titanium  exhibited  signs  of  adiabatic  shear  banding  in  associated  U.S. 
Army  Research  Laboratory  Survivability/Lethality  Analysis  Directorate  experiments,  and  this 
damage  mechanism  was  not  represented  in  the  Johnson-Cook  fracture  model. 

If  the  steel  residual  hits  the  bed  of  propellant  with  a  substantial  velocity,  effects  occur  that  were 
neglected  in  this  analysis.  The  work  performed  on  the  propellant  would  give  rise  to  additional 
heating  that  could  contribute  to  combustion  initiation. 
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Appendix.  Heat  Condition  Between  Two  Semi-Infinite  Layers,  Each  Initially 
at  a  Uniform  Layer 


A.l  Problem  Statement 

Consider  two  homogeneous,  isotropic  materials  that  obey  Fourier’s  law.  Material  1  is 
characterized  by  density  p\,  constant-pressure  specific  heat  cpi,  and  thermal  conductivity  K\. 
Material  2  is  characterized  by  pi,  cp 2,  and  i<2-  These  conductivities,  densities,  and  specific  heats 
are  assumed  to  be  material  constants,  independent  of  temperature.  The  thermal  diffusivities  of 
the  two  materials,  a\  and  «2,  are  defined  by 


ax 


P\  cpi 


(A-l) 


and 


Initially,  material  1  is  at  uniform  temperature  Qo  and  material  2  is  at  zero  temperature.  The  initial 
conditions  are  then 

6X (x,0)  =  0O ;  x<0  (A-5) 

and 

02  (x,0)  =  0 ;  x  >  0 .  (A-6) 
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Figure  A- 1 .  Problem  definition. 


Along  the  interface  between  steel  and  propellant,  perfect  thermal  contact  is  assumed  so  that 
temperature  and  heat  flux  are  both  continuous.  The  interface  boundary  conditions  are 

0{(O,t)  =  02(O,t);  t>0  (A-7) 

and 


Since  the  domains  of  materials  1  and  2  are  semi-infinite,  the  remaining  boundary  conditions  are 


0X  (x,  t)  — ^  0q  as  x  — ^  — go 


(A-9) 


and 


A.2  Analytical  Solution 


0 2  (x,  t)  — ^  0  as  x  — ^  +oo . 


(A- 10) 


The  problem  can  be  solved  with  Laplace  transforms.  Let  0X  (x,  s)  and  02  (x,  5)  be  the  Laplace 
transforms  of  0X  (x,  t)  and  02  (x,  t) ,  respectively.  These  transforms  are  defined  by 
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(A-l  1) 


/v  r  oo 

0x(x,s)=\  0x(x,t)e  dt 

J  0 

and 

/•  oo 

02 (x, s)  =  02(x,t)e~s‘dt . 

J  0 


The  Laplace  transformations  of  equations  A-3  and  A-4  leads  to  the  ODEs 


d20x 

dx7 


-±0t=-^L-x<O 
ax  ax 


and 


d202 

~dx2' 


— —  6 2  =  0 ;  x  >  0 . 
a2 


(A- 12) 


(A- 13) 


(A- 14) 


The  initial  conditions,  equations  A-5  and  A-6,  have  been  absorbed  in  equations  A- 13  and  A- 14. 
Laplace  transformation  of  the  boundary  conditions  leads  to 

^(0,5)  =  #2(0,,s) ;  t  >  0,  (A-15) 


K, 


dOx  (x,  s) 


dx 


=  K~ 


d02(x,  s ) 


x  =  0 


dx 


;  t> o, 


x  =  0 


0, 


Q \  (x,  s')  — >  —  as  x  — >  —oo  3 


(A- 16) 

(A- 17) 


and 


02  (x,  t)  — ^  0  as  x  — ^  +oo . 


The  solution  of  equations  A- 13  and  A- 14  subject  to  the  boundary  conditions  is 

7  nr  A 


0x(x,s)  =  — 
s 


1 


K 


1  +  K 


V 


and 


02(x,s)  =  V“2 

l  +  K  s 


(A- 18) 


(A- 19) 


(A-20) 


where 
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The  transforms  in  equations  A- 19  and  A-20  are  inverted  using  table  F-l  in  Poulikakos.1 


KO 

0l(x,t)  =  0l o-- — 77-erfc| 
1  +  A 


r  \ 

-x 

2  yfaj 


;  x  <0 ,  t  >  0. 


Q 

01  (x,  t)  =  — —  •  erfcl 
2  1 +  K  1 


Here,  the  error  function  is  defined  by 


'  x  ^ 
2  yfa^t 


erf  (r/)  =  -^=  V e  ^ d% 

V  ’  V^Jo 


;  x  >  0 ,  t  >  0. 


(A-22) 


(A-23) 


(A-24) 


and  the  complementary  error  function  by 

erfc(// )  =  l-erf(^).  (A-25) 

Equations  A-22  and  A-23  were  evaluated  for  materials  1  and  2  associated  with  4340  steel  and 
M30A1  propellant,  respectively.  The  material  constants  were  assigned  the  values  in  table  2.  The 
results  are  shown  in  figure  A-2.  Here,  temperature  is  plotted  as  a  function  of  time  at  x  locations 
separated  by  0.625  mm.  Note  that  at  any  given  time  temperature  is  much  more  uniform  in  the 
steel  than  in  the  propellant  due  to  the  steel’s  greater  thermal  diffusivity. 

A.3  LS-DYNA  Simulation 

The  problem  was  simulated  using  LS-DYNA.  A  single  row  of  8-node  hexagonal  brick  elements 
was  used  (figure  A-3).  Each  element  was  a  cube  with  a  0.635-mm  edge  length.  These  elements 
were  geometrically  identical  to  those  used  in  the  3D  analysis  of  the  steel  residual  and  the 
propellant  grain  in  section  4.  Also,  the  same  LS-DYNA  thermal  material  model  as  in  section  4 
was  employed.  The  two  semi-infinite  domains  were  approximated  with  100  elements  each,  for  a 
total  dimension  of  63.5  mm  each. 

Figures  A-4  and  A-5  compare  the  LS-DYNA  results  with  the  analytical  solution  in  the  steel  and 
propellant,  respectively.  In  the  steel,  the  initial  cooling  is  not  as  rapid  in  the  LS-DYNA  solution 
as  in  the  analytical  solution.  Note,  however,  the  fine  temperature  scale  in  figure  A-4.  The 
solution  in  the  propellant  is  generally  more  important  in  a  computational  scheme  to  predict 
propellant  initiation,  and  the  LS-DYNA  solution  in  figure  A-5  appears  to  be  quite  good. 


^Poulikakos,  D.  Conduction  Heat  Transfer,  Prentice  Hall:  Englewood  Cliffs,  NJ,  1994. 
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Analytical  Solution 


Figure  A-2.  Analytical  results  for  temperature  vs.  time  at  the  specific  locations. 
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Figure  A-5.  LS-DYNA  results  compared  with  the  analytical  solution  in  the  M30A1  propellant. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


A,  B,  C,  M,  N 
Co, C5 
Cs 

D\,  D5 
L 
P 
S 

Us 

cp 

Cpl 

Cp2 

n\ 

n-2 

t 

Aower 

Supper 

Up 

x,  y,z 

x 

a 

a\ 

«2 

s 

s 


constants  in  the  Johnson-Cook  strength  model 
constants  in  the  Zerilli- Armstrong  FCC  strength  model 
sound  speed 

constants  in  the  Johnson-Cook  fracture  model 

rod  length 

pressure 

proportionality  constant  in  the  linear  Us-up  shock  Hugoniot  relation 
shock  speed 

specific  heat  at  constant  pressure 

constant-pressure  specific  heat  in  steel  residual 

constant-pressure  specific  heat  in  M30A1  grain 

coordinate  along  the  outward  normal  to  the  steel  residual’s  boundary 

coordinate  along  the  outward  normal  to  the  propellant  grain’s  boundary 

time 

time  below  which  sustained  combustion  does  not  occur 
time  above  which  one  can  always  expect  sustained  combustion 
particle  speed 

spatial  (laboratory)  coordinates 
position  vector 
thermal  diffusivity 
thermal  diffusivity  of  steel  residual 
thermal  diffusivity  of  M30A1  grain 
equivalent  plastic  strain 

dimensionless  time  rate  of  change  of  equivalent  plastic  strain 
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K 


*1 


K2 

V 

e 

et 


6x 

02 

dm 

0r 

P 

P\ 

Pi 

^eq 


ij  component  of  the  deviatoric  plastic  strain  tensor 

thermal  conductivity 

thermal  conductivity  of  steel  residual 

thermal  conductivity  of  M30A1  grain 

Poisson  ratio 

temperature 

homologous  temperature 

CTH-determined  temperature  field  in  steel  residual 

temperature  field  in  steel  residual 

temperature  field  in  M30A1  grain 

melt  temperature 

room  temperature 

density 

density  of  steel  residual 

density  of  M30A1  grain 

von  Mises  Cauchy  stress 

ij  component  of  the  Cauchy  stress  tensor 

ij  component  of  the  deviatoric  Cauchy  stress  tensor 
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